Fluctuation broadening of plasma resonance line in the vortex liquid state of layered 

superconductors 
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The Josephson plasma resonance (JPR) provides a sensitive probe of vortex states in layered 
superconductors. We demonstrate that in the case of weak damping in the liquid phase, broadening 
of the JPR line is caused mainly by random Josephson coupling arising from the density fluctuations 
of pancake vortices. In this case the JPR line has the universal shape, which is determined only by 
parameters of the superconductors and temperature. This mechanism gives a natural explanation 
for the experimentally observed asymmetric lineshape. The tail at high frequencies arises due to 
mixing of the propagating plasma modes by random Josephson coupling, while the tail at small 
frequencies is caused by the localized plasma modes originating from a rare fluctuation suppression 
of the Josephson coupling in large areas. 



Observation of the Josephson plasma resonance (JPR) 
in the mixed state of iawered high temperature and or- 
ganic superconductoraJ □ allows the study of Josephson 
coupling of layers over a wide range of fields and tem- 
peratures. The Josephson coupling in the mixed state is 
mainly determined by correlations between arrangements 
of pancake vortices in neighboring layers and therefore 
it can serve as a probe of vortex state. The resonance 
absorption is determined by the spatial and temporal 
behavior of the "local coherence parameter" C„(r,t) = 
cos(^„,„+i(r,<), with (p„^„+i = (pn+i - (fin - (27rs/$o)^2 
being the gauge-invariant phase difference between the 
layers n and n + 1. Here r = {x,y) and z are the 
coordinates in the ab plane and along the c axis, s is 
the interlayer spacing. Nonzero ipn^n+i arises due to the 
misalignment of pancake vortices in neighboring layers. 
Time variations of ipn,n+i caused by pancake motion are 
usually much slower than the plasma oscillations. In 
this case the plasma mode probes a snapshot of the in- 
stantaneous phase distribution. In the liquid phase at 
high fields, B ^ Bj, correlations between pancake po- 
sitions in neighboring layers are almost absent. Here 
Bj ~ <i>o/Aj, Aj — js is the Josephson length, and 7 
is the anisotropy factor. In this state C„(r) rapidly os- 
cillates in space so that its variations are much larger 
than average. The resonance in such a situation occurs 
because small phase oscillations induced by the exter- 
nal electric field change slowly in space and average out 
these rapid variations. We will estimate that the reso- 
nance is formed at the typical scale Aj/a, with a being 
the intervortex spacing. This large scale averaging leads 
to a fairly sharp resonance, with the resonance frequency 
squared ujp being roughly proportional to the averaged 
cosine factor C = C(B,T) (see, e.g., Refs. |,|) 



cul{B,T) « ujI{T)C, C = (cos^„,„+i(r)) 



(1) 



and (. . .) denotes the thermodynamic average. However, 
the averaging by the smooth oscillating phase is not com- 
plete. Large scale fluctuations of C„(r), induced by pan- 
cake density fluctuations, lead to inhomogeneous broad- 
ening of the JPR line. This mechanism was proposed 
in Ref. ^ to explain line broadening in the vortex glass 
state. In this paper we analyze broadening of the JPR 
line in the vortex liquid due to the random Josephson 
coupling. An attractive feature of the line broadening in 
the liquid phase is that it is an intrinsic property of the 
material caused by thermal fluctuations. The plasma res- 
onance line has a very peculiar asymmetric shape with 
the long tail in the high field partErEI. This shape can 
be naturally explained by the proposed mechanism. The 
high frequency/high field tail of the line probes mixing of 
the propagating plasma modes by the random Josephson 
coupling, while the low frequency/low field tail probes 
the localized plasma modes originating due to rare fluc- 
tuation suppression of the coupling. We will analytically 
calculate the resonant absorption in both regions. 

The external alternating electric field applied along the 
c axis Vz — 'Dzo exp(ja;t) induces small oscillations of the 
interlayer phase difference ^'n^n+i- general, charging 
of layersQ and deviations o£,thc quasiparticle distribution 
function from equilibriumu should be accounted for in 
the time-dependent equation for the phase difference. To 
illustrate the physics of the inhomogeneous line broaden- 
ing, we will consider a simplified, equation, in which we 
do not account for these effect^Ha 



-2 + ^'jLV' 



iuj'Dz 



(2) 



where Lo^iT) — /ecX^iT), Ac is the c-component of the 
London penetration depth, Cc is the dielectric constant. 



where Jq — c<I>o/87r^A^s is the Josephson current. The 
second term in the brackets accounts for the induc- 
tive interaction between the junctions. The operator 
L acts on the layer index n and is defined as LAn ~ 
J2„iLnniAn with L„„, = J^"" ^ L{q) cos{n - m)q, and 
L{q) = [2(1 — cos (7) + s^/A^f,] , with Xat being the ab- 
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components of the London penetration depth. We ne- 
glect in Eq. (j^) time variations of C„(r, t) assuming them 
to be small during the time l/ut. If we neglect the charg- 
ing effectaZl, then the oscillating internal electric field is 
connected with the solution of Eq. (H) by the Josephson 
relation Ez ~ (ia;$o/27rcs)(p'j The resonant absorp- 
tion is given by the imaginary part of the inverse dielec- 
tric constant l/£c(w) = Ez/Vz- We split C„(r) into the 
average and fluctuating parts, C„(r) = C + Un{r). The 
correlation function of u„(r) is given by 

{un{r)un'{r')) = (cos(^„^„+i(r)cos.^„'^„'+i(r')) - 

This correlation function depends weakly on Ej. Taking 
the limit Ej — > 0, we obtain 

(M„(r)u„/(r')) sa ^S{r - r')5nn' , 

where S{r) — (cos[iy9„.„+i(r) — </?„,„+i(0)]) is the static 
phase correlation function at Ej = 0. Static configu- 
rations ipn,n+i(j) are mainly determined by the thermal 
fluctuations of pancake vortices, and S{r) drops at dis- 
tances of the order of the intervortex spacing a. 

An important observation is that in spite of the 
rapid variations of C„(r) in Eq. (||), the solution ip' = 
^p'n n+ii''^) varies smoothly in space. At u ^ Up the 
typical length scale of phase variations can be es- 
timated by balancing the typical kinetic energy of su- 
percurrents, Eo{(p')'^ / L"^, with the typical fluctuation 
of the random Josephson energy, Ej{(p')'^a/L^. Here 
Eo = s<f>Q/167r^A^j is the in-plane phase stiffness and 
Ej = Eq/X^j is the Josephson energy per unit area. 
This gives — Xj/a 3> a. Smoothly varying ip' ef- 
fectively averages the rapid variation of C„(r) and the 
plasma frequency is simply determined by C(B, T). Fluc- 
tuations of C„(r) smoothened over the area L^, Cl^ — 
^ip^ Ir<L '^^^ni'i'^), produce inhomogeneous broadening 
of the JPR line. Calculating the mean squared fluctua- 
tion of Ci^, {{Cl, - Cf) « aVL| = a^/A^, we obtain 
the estimate for the inhomogeneous line width 



■Lol{T)B.,/B. 



(3) 



Inhomogeneous line broadening is determined by the 
amplitude of fluctuations of w„(r), Sq — J drS{r) w a^, 
i.e., cx 5*0. On the OtJier hand, according to the high 
temperature expansionoEj, the average cosine is also de- 
termined by Sq: 



CiB,T)^{Ej/2T)So, 



(4) 



This relation allows us to connect the linewidth with the 
resonance frequency. In Ref. |o| 5*0 and C(B, T) were cal- 
culated taking into account both the vortex and regular 
phase fluctuations. At low temperatures, T <C ^ttEq, the 
result can be written as 



C(B,T)«M2|£,i 



2^i;o \Bj 



(5) 



where fs{T) is the universal function ofT /Eq. The loga- 
rithmic correction appears due to the regular phase fluc- 
tuations. Comparing Eqs. (Q) and (||) with Eq. (||), 
we obtain an estimate for the relative line width in 
the liquid state due to the inhomogeneous broadening, 

LOb/uop kT/Eq <t: 1. 

Consider now the problem quantatively. Eq. (|^) is sim- 
ilar to the Schrodinger equation for a particle in a random 
potential and the problem is mapped onto the problem of 
the density of states in such a system. Then one can use 
the same techniques. The problem of density of states 
does not have an analytical solution over the whole en- 
ergy range. Only asymptotics for positive and negative 
energies can be calculated. 

At frequencies well above the resonance frequency one 
can treat the random coupling Un(r) as a perturbation. 
A perturbative analyses of Eq.(^ can be conveniently 
performed using the Green function formalism. We de- 
fine the Green function Qnn' (r, r'; E) as a solution of the 
equation. 



E + X^jLV^ - unir) Gnn' (r, r'; E) = S^n'S (r - r') 



Here we introduced the dimensionless "energy" E = 
uj'^/ujq — C. Knowledge of Gnn' (r, r'; E) allows one to cal- 
culate the phase responses to arbitrary external pertur- 
bations. In particular, the response to the homogeneous 
external field is given by the averaged Fourier transform 
oiGnn' {r,r';E), 



g (k, q;E) J '^'^ (^"0 (r, 0; E)) exp {~ikr - 



iqn) , 



at k, (7 = 0. The dielectric constant ec{(^) is connected 
with g (k, q; E) as ejedto) = {u? Idl)g (O; wVwg - C) . 
The JPR line, given by the imaginary part of 1 / ec{'jj), can 
be connected with the spectral density A[){E), Aq{E) = 
Im[(5 {{)-E-i5))]/T: 



p{uj) = Im 



LU, 







LU, 







(6) 



The perturbative expansion of g (k, q; E) with respect 
to u„{r) can be performed using the diagrammatic tech- 
nique and g (k, q; E) is represented as 
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E - L{q)Vjk'^ -ll{\<i,q-E)' 



(7) 



where S (k, q] E) is the self energy function. It, in turn, 
can be represented as the perturbation series with respect 
to u„(r). In the lowest order with respect to w„(r), which 
is equivalent to the Born approximation for scattering, 
S (k, q\ E), is given by 



^{k,q;E) 



^C;(k',g';i?)5(k-k') (8) 



The imaginary part of E(k, (7;i?), which we denote by 
S2, determines the line broadening, while the real part. 
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Si, determines the shift of the resonance frequency due 
to inhomogeneities. The high frequency asymptotics 
of E (k, q\ E) can be obtained by replacing Q (k, q\ E) 



with its bare value Qq (k, q; E) = \E — L{q)X'jk'^ 



isy 



The main contribution to S2 comes from the region 
k ~ \/E/Xj ^ 1/a. In this region one can neglect the 
/c-dependence of S(k) and replace it by Sq. The same 
replacement is possible for almost all terms in the per- 
turbation series. The only exception is the Born integral 
for Si logarithmically diverging at fc ^ ^/E/\j. This di- 
vergency is cut by the k dependence of ^(k) at /c « 1/a. 
Using the replacement S'(k) — > 5o ~ {2T/Ej)C and per- 
forming integrations with respect to k and q we obtain 
the following expressions for the imaginary and real parts 
of the self energy function at large "energies" 



-'2\E- 



TC 
2^' 



Sil 



m „ ■ 

TT a^E 



(9) 



Soo ~ l^j provides the typical "energy" scale for 
the broadening. In Si we separate the logarithmic 
contribution, — (2Soo/7r) ln(Aj/a^), coming from large 
k, and the remaining part which is determined by Sq. 
Further analysis of the perturbation expansion shows 
that the expansion parameter at large E is Soo/i-B + 
(2Soo/7r) ln(Aj/a^)]. Therefore the components of the 
self energy can be represented in the following scaling 
form 



Si 



A2 



^=S2(C), -^ = — In^-.i(C) 



(10) 



with C = -E/Soo + (2/7r)ln(Aya2). si(C) and s-^iC,) are 
the dimensionless functions with the asymptotics, S2 — > 1 
and S\ (I/tt) ln(l/^) at ^ 00. Using this expressions 
we can represent the JPR absorption dq) in scaling form 

p(a;KVc^2 = /(C), c = (^'-<:^p')M' (11) 



/(C) = 



«2(C) 



[C + si(C)]^ + -'^i(C)' 



(12) 



where Wh is the resonance width, lj^ — {T/2Eo)ujp, 
in agreement with the estimate (||). We assume this 
width to be larger than the linewidth due to the quasi- 



particle dissipation. 



is the resonance frequency 



shifted by the random Josephson coupling, uj^ = u!p{l — 
{T/TTEo)ln{B/Bj)). This negative shift is due to the 
second order perturbative correction to the ground state 
"energy" , similar to the well known result of quantum 
mechanical perturbation theory. Note that the inho- 
mogeneous correction is two times larger than the cor- 
rection due to the regular phase fluctuations (see Eq. 
(^) and has the opposite sign. The high frequency tail 
of the resonant absorption, oj^ — aj^ ^ w^, is given 

by p{ll>) ~ ^'^^b/ (^^ ~'^p)^- The whole shape of the 
line is determined by the single dimensionless parameter 
/X = T/ Eq. In particular, the width of the line is /xo-ip, the 
maximum absorption scales as and the absorption 



at Lu"^ — Ki scales as ^. Using the field dependence 
of the plasma frequency given by Eqs. (|l|) and (|^) we can 
also present the scaling parameter C as a function of the 
magnetic field 



C=iB- Br) /Bb + (I/tt) HB/Bj), 



(13) 



with Br = /.cj2^.7$o/(2Tu;2) ^nd Bb ^ ifsU;^/Au^)Bj^ 
{T/2Eo)Br. This allows us to obtain from Eqs. (0) 
and (O) the field dependence of the resonant absorption, 
which is usually probed in the JPR experiments. 

The problem is now reduced to a calculation of the 
dimensionless functions si and S2. These functions can 
be calculated approximately if wc keep only the first term 
in the expansion of S given by Eq. (||) (self consistent 
Born approximation (SCBA)). This leads to the following 
equations 



11 , C + ^i 

S2 = — ! — arctan 

2 TT S2 

si = -^ln[(C + si)2 + . 



(14) 



which we solve numerically. Fig. |l] shows the dependence 
of the reduced JPR absorption p(ti;)(jj^ /oj^ on the reduced 
frequency C- Note that the SCBA actually breaks down 
at ^ ^ 1, i.e., it describes quantitatively only the right 
hand side of the line. For comparison we took the exper- 
imental JPR line for optimally dopjsd Bi2Sr2CaCu208 
(BSCCO) crystal with Tc = 89.5KEI at u;/27r = 45GHz 
and T — 51. 4K. We replotted this line as the function of 
the scaling parameter ^ from Eq. (p^), where the param- 
eters Br and Bb are chosen to make the experimental and 
theoretical lines match at large C giving Br = 0.053T and 
Bb = 0.016T. The obtained ratio Bb/Br ~ 0.3 is some- 
what larger then the theoretical estimate T/2Eq « 0.17, 
which wc obtain taking Xab = 200nm/y^l - (T/Tc)'^. 

The SCBA approach gives the JPR line terminating 
at some finite frequency. In reality a long absorption tail 
exists on the low frequency side of the line due to the 
rare fiuctuation configurations corresponding to suppres- 
sion of C„(r) over large areas. This is very similar to 
the localization tail below the bottom of the conduction 
band in disordered semiconductors (Lifshitz tail, see, e.g., 
Ref. |l^). An important simplification is that due to the 
large scale of spatial changes of </'^„+i(?'), the random 
"potential" u„ (r) can be treated as a short ranged Gaus- 
sian random variable with the probability distribution 



P[u„(r)] cx exp [-n{Un}] , n{Un} dv 



K(r)]' 
5*0 



Following a standard line of reasoning, we estimate 
the spectral density Aq[E) with exponential accuracy as 
Aq{E) (X exp(-$(i;)) with 



(15) 



where £q{u\ is the ground state energy for a given po- 
tential fluctuation M„(r), 
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£o{u} — min 



i/o{*} 



E 



dru„(r)*2(r) 



(16) 



with ffo{*} = ^„ „^ / drL„„V*„V*,„ and normal- 
ization ^Ylin I '^^^n ~ ^- ^ conditional minimization in 
Eq. ( p^ can be performed using the Lagrange technique 



= -pE- 



(17) 



mm 



n{un{r)} 



where /3 is the Lagrange factor, which has to be found 
from the relation £o{u} = E. Optimization with respect 
to ^^(r) gives u„(r) = — S'o/3^'^(r)/2. Substituting this 



expression into Eq. (17) and variating it with respect to 
we obtain the nonlinear eigenvalue problem 

- \^jLVHr.{v) - (/35o/2)*3(r) = i?*„(r), (18) 

which determines the optimum and (3. To sim- 
plify this equation we (i) apply the operator L^^ — 
s^/A^^ + both sides and neglect small term s^/A^^ 

in the resulting equation (here V^j is defined as V^Q!„ = 
Q!„+i — 2q!„) and (ii) introduce dimensionless vari- 
ables Vn(r) = «'„(r)/*o(0), f = rV/55'o^'^(0)/2/Aj, and 
to = 2£'/(/3S'o\E'o(0)). These transformations lead to the 
dimensionless nonlinear eigenvalue problem 

-eoVliJnir), 



(19) 



which has to be solved with the condition V'o (0) = 1 and 
with the dipolelike asymptotics at large distances 

i^„(r) = -aV^deol + n^)~^'\ f, n » 1. 

The spectral density Ao{E) can be expressed through 
and V'ii(f) 

Ao{E) (X exp (-^||^) , ^4 = E / dii'n- (20) 

A numerical solution of Eq. ( p^ gives to — —0.164, ip^ = 
2.418, and a = 0.224. From Eq. dH) we now obtain 
Ao{E) cx exp (— 14.72Aji?/5o), which corresponds to the 
exponential tail of the JPR absorption 



p{lli) cx exp(3.18C) = exp 



7.36£:o {^p - 



(21) 



An exponential fit of the left hand side of the experimen- 
tal line in Fig. |l] gives p oc exp(3.63C), consistent with 
above result. 

In conclusion, we studied the universal shape of the 
JPR line due to thermal fluctuations of pancakes in the 
vortex liquid. We showed that the line broadening is 
strongly asymmetric, in agreement with experimental 
observations. We found that the relative linewidth is 



given by the parameter T/Eq, which also determines the 
strength of thermal fluctuations at i? = and is directly 
related to the temperature dependent Xab- 
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FIG. 1. 



Plot of the reduced JPR lineshape p{LL})ujh I'-o'^ 



reduced frequency C, obtained within SCBA (Eqs. plh and 
(^^). For comparison we also show the experimental JPR 
line for BSCCO plotted as the function of C, from Eq. ( |l^ ) 
(courtesy of M. GaifuUin and Y. Matsuda, inset shows raw 
data). The parameters Br and Bi, in Eq. (jl^) are chosen to 
make the experimental and theoretical lines match at large C, 
giving Br = 0.053T and B^ = 0.016T. The left hand side of 
the line is fitted to exponent for comparison with Eq. d2ll) 
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